Alternative splicing analysis of lignocellulose-degrading enzyme genes and enzyme variants in Aspergillus niger.

Alternative splicing (AS) greatly expands the protein diversity in eukaryotes. Although AS variants have been frequently reported existing in filamentous fungi, it remains unclear whether lignocellulose-degrading enzyme genes in industrially important fungi undergo AS events. In this work, AS events of lignocellulose-degrading enzymes genes in Aspergillus niger under two carbon sources (glucose and wheat straw) were investigated by RNA-Seq. The results showed that a total of 23 out of the 56 lignocellulose-degrading enzyme genes had AS events and intron retention was the main type of these AS events. The AS variant enzymes from the annotated endo-ß-1,4-xylanase F1 gene (xynF1) and the endo-ß-1,4-glucanase D gene (eglD), noted as XYNF1-AS and EGLD-AS, were characterized compared to their normal splicing products XYNF1 and EGLD, respectively. The AS variant XYNF1-AS displayed xylanase activity whereas XYNF1 did not. As for EGLD-AS and EGLD, neither of them showed annotated endo-ß-1,4-glucanase activity. Instead, both showed lytic polysaccharide monooxygenase (LPMO) activity with some differences in catalytic properties. Our work demonstrated that the AS variants in A. niger were good sources for discovering novel lignocellulose-degrading enzymes. KEY POINTS: • AS events were identified in the lignocellulose-degrading enzyme genes of A. niger. • New ß-1,4-xylanase and LPMO derived from AS events were characterized.

PRPF8-mediated dysregulation of hBrr2 helicase disrupts human spliceosome kinetics and 5´-splice-site selection causing tissue-specific defects.

The carboxy-terminus of the spliceosomal protein PRPF8, which regulates the RNA helicase Brr2, is a hotspot for mutations causing retinitis pigmentosa-type 13, with unclear role in human splicing and tissue-specificity mechanism. We used patient induced pluripotent stem cells-derived cells, carrying the heterozygous PRPF8 c.6926 A > C (p.H2309P) mutation to demonstrate retinal-specific endophenotypes comprising photoreceptor loss, apical-basal polarity and ciliary defects. Comprehensive molecular, transcriptomic, and proteomic analyses revealed a role of the PRPF8/Brr2 regulation in 5'-splice site (5'SS) selection by spliceosomes, for which disruption impaired alternative splicing and weak/suboptimal 5'SS selection, and enhanced cryptic splicing, predominantly in ciliary and retinal-specific transcripts. Altered splicing efficiency, nuclear speckles organisation, and PRPF8 interaction with U6 snRNA, caused accumulation of active spliceosomes and poly(A)+ mRNAs in unique splicing clusters located at the nuclear periphery of photoreceptors. Collectively these elucidate the role of PRPF8/Brr2 regulatory mechanisms in splicing and the molecular basis of retinal disease, informing therapeutic approaches.

Loss-of-function mutation in PRMT9 causes abnormal synapse development by dysregulation of RNA alternative splicing.

Protein arginine methyltransferase 9 (PRMT9) is a recently identified member of the PRMT family, yet its biological function remains largely unknown. Here, by characterizing an intellectual disability associated PRMT9 mutation (G189R) and establishing a Prmt9 conditional knockout (cKO) mouse model, we uncover an important function of PRMT9 in neuronal development. The G189R mutation abolishes PRMT9 methyltransferase activity and reduces its protein stability. Knockout of Prmt9 in hippocampal neurons causes alternative splicing of ~1900 genes, which likely accounts for the aberrant synapse development and impaired learning and memory in the Prmt9 cKO mice. Mechanistically, we discover a methylation-sensitive protein-RNA interaction between the arginine 508 (R508) of the splicing factor 3B subunit 2 (SF3B2), the site that is exclusively methylated by PRMT9, and the pre-mRNA anchoring site, a cis-regulatory element that is critical for RNA splicing. Additionally, using human and mouse cell lines, as well as an SF3B2 arginine methylation-deficient mouse model, we provide strong evidence that SF3B2 is the primary methylation substrate of PRMT9, thus highlighting the conserved function of the PRMT9/SF3B2 axis in regulating pre-mRNA splicing.

SF3A2 promotes progression and cisplatin resistance in triple-negative breast cancer via alternative splicing of MKRN1.

Triple-negative breast cancer (TNBC) is the deadliest subtype of breast cancer owing to the lack of effective therapeutic targets. Splicing factor 3a subunit 2 (SF3A2), a poorly defined splicing factor, was notably elevated in TNBC tissues and promoted TNBC progression, as confirmed by cell proliferation, colony formation, transwell migration, and invasion assays. Mechanistic investigations revealed that E3 ubiquitin-protein ligase UBR5 promoted the ubiquitination-dependent degradation of SF3A2, which in turn regulated UBR5, thus forming a feedback loop to balance these two oncoproteins. Moreover, SF3A2 accelerated TNBC progression by, at least in part, specifically regulating the alternative splicing of makorin ring finger protein 1 (MKRN1) and promoting the expression of the dominant and oncogenic isoform, MKRN1-T1. Furthermore, SF3A2 participated in the regulation of both extrinsic and intrinsic apoptosis, leading to cisplatin resistance in TNBC cells. Collectively, these findings reveal a previously unknown role of SF3A2 in TNBC progression and cisplatin resistance, highlighting SF3A2 as a potential therapeutic target for patients with TNBC.

Comprehensive deciphering the alternative splicing patterns involved in leaf morphogenesis of Liriodendron chinense.

Alternative splicing (AS), a pivotal post-transcriptional regulatory mechanism, profoundly amplifies diversity and complexity of transcriptome and proteome. Liriodendron chinense (Hemsl.) Sarg., an excellent ornamental tree species renowned for its distinctive leaf shape, which resembles the mandarin jacket. Despite the documented potential genes related to leaf development of L. chinense, the underlying post-transcriptional regulatory mechanisms remain veiled. Here, we conducted a comprehensive analysis of the transcriptome to clarify the genome-wide landscape of the AS pattern and the spectrum of spliced isoforms during leaf developmental stages in L. chinense. Our investigation unveiled 50,259 AS events, involving 10,685 genes (32.9%), with intron retention as the most prevalent events. Notably, the initial stage of leaf development witnessed the detection of 804 differentially AS events affiliated with 548 genes. Although both differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) were enriched into morphogenetic related pathways during the transition from fishhook (P2) to lobed (P7) leaves, there was only a modest degree of overlap between DASGs and DEGs. Furthermore, we conducted a comprehensively AS analysis on homologous genes involved in leaf morphogenesis, and most of which are subject to post-transcriptional regulation of AS. Among them, the AINTEGUMENTA-LIKE transcript factor LcAIL5 was characterization in detailed, which experiences skipping exon (SE), and two transcripts displayed disparate expression patterns across multiple stages. Overall, these findings yield a comprehensive understanding of leaf development regulation via AS, offering a novel perspective for further deciphering the mechanism of plant leaf morphogenesis.

The landscape of abiotic and biotic stress-responsive splice variants with deep RNA-seq datasets in hot pepper.

Alternative splicing (AS) is a widely observed phenomenon in eukaryotes that plays a critical role in development and stress responses. In plants, the large number of RNA-seq datasets in response to different environmental stressors can provide clues for identification of condition-specific and/or common AS variants for preferred agronomic traits. We report RNA-seq datasets (350.7 Gb) from Capsicum annuum inoculated with one of three bacteria, one virus, or one oomycete and obtained additional existing transcriptome datasets. In this study, we investigated the landscape of AS in response to environmental stressors, signaling molecules, and tissues from 425 total samples comprising 841.49 Gb. In addition, we identified genes that undergo AS under specific and shared stress conditions to obtain potential genes that may be involved in enhancing tolerance to stressors. We uncovered 1,642,007 AS events and identified 4,354 differential alternative splicing genes related to environmental stressors, tissues, and signaling molecules. This information and approach provide useful data for basic-research focused on enhancing tolerance to environmental stressors in hot pepper or establishing breeding programs.

Genome-wide characterization of post-transcriptional processes related to wood formation in Dalbergia odorifera.

BACKGROUND: Alternative polyadenylation (APA), alternative splicing (AS), and long non-coding RNAs (lncRNAs) play regulatory roles in post-transcriptional processes in plants. However, little is known about their involvement in xylem development in Dalbergia odorifera, a valuable rosewood species with medicinal and commercial significance. We addressed this by conducting Isoform Sequencing (Iso-Seq) using PacBio's SMRT technology and combined it with RNA-seq analysis (RNA sequencing on Illumina platform) after collecting xylem samples from the transition zone and the sapwood of D. odorifera. RESULTS: We identified 14,938 full-length transcripts, including 9,830 novel isoforms, which has updated the D. odorifera genome annotation. Our analysis has revealed that 4,164 genes undergo APA, whereas 3,084 genes encounter AS. We have also annotated 118 lncRNAs. Furthermore, RNA-seq analysis identified 170 differential alternative splicing (DAS) events, 344 genes with differential APA site usage (DE-APA), and 6 differentially expressed lncRNAs in the transition zone when compared to the sapwood. AS, APA, and lncRNAs are differentially regulated during xylem development. Differentially expressed APA genes were enriched for terpenoid and flavonoid metabolism, indicating their role in the heartwood formation. Additionally, DE-APA genes were associated with cell wall biosynthesis and terpenoid metabolism, implying an APA's role in wood formation. A DAS gene (involved in chalcone accumulation) with a significantly greater inclusion of the last exon in the transition zone than in the sapwood was identified. We also found that differentially expressed lncRNAs targeted the genes related to terpene synthesis. CONCLUSIONS: This study enhances our understanding of the molecular regulatory mechanisms underlying wood formation in D. odorifera, and provides valuable genetic resources and insights for its molecular-assisted breeding.

High-quality faba bean reference transcripts generated using PacBio and Illumina RNA-seq data.

The genome of faba bean was first published in 2023. To promote future molecular breeding studies, we improved the quality of the faba genome based on high-density genetic maps and the Illumina and Pacbio RNA-seq datasets. Two high-density genetic maps were used to conduct the scaffold ordering and orientation of faba bean, culminating in an increased length (i.e., 14.28 Mbp) of chromosomes and a decrease in the number of scaffolds by 45. In gene model mining and optimisation, the PacBio and Illumina RNA-seq datasets from 37 samples allowed for the identification and correction 121,606 transcripts, and the data facilitated a prediction of 15,640 alternative splicing events, 2,148 lncRNAs, and 1,752 fusion transcripts, thus allowing for a clearer understanding of the gene structures underlying the faba genome. Moreover, a total of 38,850 new genes including 56,188 transcripts were identified compared with the reference genome. Finally, the genetic data of the reference genome was integrated and a comprehensive and complete faba bean transcriptome sequence of 103,267 transcripts derived from 54,753 uni-genes was formed.

Purification of SRSF1 from E. coli for Biophysical and Biochemical Assays.

The Ser/Arg-rich splicing factors (SR proteins) constitute a crucial protein family in alternative splicing, comprising twelve members characterized by unique repetitive Arg-Ser dipeptide sequences (RS) and one to two RNA-recognition motifs (RRM). The RS regions of SR proteins undergo variable phosphorylation, resulting in unphosphorylated, partially phosphorylated, or hyper-phosphorylated states based on functional requirements. Despite the identification of the SR protein family over 30 years ago, the purification of native SR proteins in soluble form at large quantities has presented challenges due to their low solubility. This protocol delineates a method for acquiring soluble, full-length, unphosphorylated, hypo- and hyper-phosphorylated SRSF1, a prototypical SR family member. Notably, this protocol facilitates the purification of SRSF1 in ample quantities suitable for NMR, as well as various biophysical and biochemical studies. The methodologies and principles outlined herein are expected to extend beyond SRSF1 protein production and can be adapted for purifying other SR protein family members or SR-related proteins, such as snRNP70 and U2AF-35. Given the involvement of these proteins in numerous essential biological processes, this protocol will prove beneficial to researchers in related fields. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Purification of SRSF1 from E. coli Support Protocol: Purification of ULP1 Basic Protocol 2: Purification of hypo-phosphorylated SRSF1 from E. coli Basic Protocol 3: Purification of hyper-phosphorylated SRSF1 from E. coli.

SRRM2 splicing factor modulates cell fate in early development.

Embryo development is an orchestrated process that relies on tight regulation of gene expression to guide cell differentiation and fate decisions. The Srrm2 splicing factor has recently been implicated in developmental disorders and diseases, but its role in early mammalian development remains unexplored. Here, we show that Srrm2 dosage is critical for maintaining embryonic stem cell pluripotency and cell identity. Srrm2 heterozygosity promotes loss of stemness, characterised by the coexistence of cells expressing naive and formative pluripotency markers, together with extensive changes in gene expression, including genes regulated by serum-response transcription factor (SRF) and differentiation-related genes. Depletion of Srrm2 by RNA interference in embryonic stem cells shows that the earliest effects of Srrm2 heterozygosity are specific alternative splicing events on a small number of genes, followed by expression changes in metabolism and differentiation-related genes. Our findings unveil molecular and cellular roles of Srrm2 in stemness and lineage commitment, shedding light on the roles of splicing regulators in early embryogenesis, developmental diseases and tumorigenesis.

Antisense oligonucleotide therapeutic approach for Timothy syndrome.

Timothy syndrome (TS) is a severe, multisystem disorder characterized by autism, epilepsy, long-QT syndrome and other neuropsychiatric conditions1. TS type 1 (TS1) is caused by a gain-of-function variant in the alternatively spliced and developmentally enriched CACNA1C exon 8A, as opposed to its counterpart exon 8. We previously uncovered several phenotypes in neurons derived from patients with TS1, including delayed channel inactivation, prolonged depolarization-induced calcium rise, impaired interneuron migration, activity-dependent dendrite retraction and an unanticipated persistent expression of exon 8A2-6. We reasoned that switching CACNA1C exon utilization from 8A to 8 would represent a potential therapeutic strategy. Here we developed antisense oligonucleotides (ASOs) to effectively decrease the inclusion of exon 8A in human cells both in vitro and, following transplantation, in vivo. We discovered that the ASO-mediated switch from exon 8A to 8 robustly rescued defects in patient-derived cortical organoids and migration in forebrain assembloids. Leveraging a transplantation platform previously developed7, we found that a single intrathecal ASO administration rescued calcium changes and in vivo dendrite retraction of patient neurons, suggesting that suppression of CACNA1C exon 8A expression is a potential treatment for TS1. Broadly, these experiments illustrate how a multilevel, in vivo and in vitro stem cell model-based approach can identify strategies to reverse disease-relevant neural pathophysiology.

FOXP3 full length splice variant is associated with kidney allograft tolerance.

Background: Progressive decline of allograft function leads to premature graft loss. Forkhead box P3 (FOXP3), a characteristic gene of T-regulatory cells, is known to be essential for auto-antigen tolerance. We assessed the hypothesis that low FOXP3 mRNA splice variant levels in peripheral blood cells early after transplantation are associated with progressive allograft injury. Methods: Blood samples were prospectively collected from 333 incident kidney transplant recipients on the first and 29th postoperative day. We used quantitative polymerase chain reaction to determine transcripts of 3 isotypes of FOXP3 splice variants, including pre-mature FOXP3 and full length FOXP3 (FOXP3fl). We investigated the association between FOXP3 splice variant levels and the declines in estimated glomerular filtration rate (eGFR) of more than 5ml/min/1.73m2 within the first-year post-transplant using logistic regression. Results: We observed lower FOXP3fl levels in recipients with declining eGFR (N = 132) than in recipients with stable eGFR (N = 201), (logarithmic value -4.13 [IQR -4.50 to -3.84] vs -4.00 [4.32 to -3.74], p=0.02). In ad hoc analysis pre-transplant FOXP3fl levels were similar in both groups. The association between FOXP3fl and declining eGFR was confirmed by multivariable analysis adjusted for potential confounding factors (Odds Ratio 0.51, 95% confidence interval 0.28 to 0.91: p=0.02). When stratifying FOXP3fl levels into quartiles, recipients with lower day1 FOXP3fl had the highest rate of declining eGFR (p=0.04). Conclusion: Low FOXP3fl splice variant levels at the first postoperative day in kidney transplant recipients were associated with severe decline of eGFR, a well-known surrogate for hard endpoints.

Alternative splicing for tuneable expression of protein subunits at desired ratios.

The controlled expression of two or more proteins at a defined and stable ratio remains a substantial challenge, particularly in the bi- and multispecific antibody field. Achieving an optimal ratio of protein subunits can facilitate the assembly of multimeric proteins with high efficiency and minimize the production of by-products. In this study, we propose a solution based on alternative splicing, enabling the expression of a tunable and predefined ratio of two distinct polypeptide chains from the same pre-mRNA under the control of a single promoter. The pre-mRNA used in this study contains two open reading frames situated on separate exons. The first exon is flanked by two copies of the chicken troponin intron 4 (cTNT-I4) and is susceptible to excision from the pre-mRNA by means of alternative splicing. This specific design enables the modulation of the splice ratio by adjusting the strength of the splice acceptor. To illustrate this approach, we developed constructs expressing varying ratios of GFP and dsRED and extended their application to multimeric proteins such as monoclonal antibodies, achieving industrially relevant expression levels (>1 g/L) in a 14-day fed-batch process. The stability of the splice ratio was confirmed by droplet digital PCR in a stable pool cultivated over a 28-day period, while product quality was assessed via intact mass analysis, demonstrating absence of product-related impurities resulting from undesired splice events. Furthermore, we showcased the versatility of the construct by expressing two subunits of a bispecific antibody of the BEAT® type, which contains three distinct subunits in total.

Unveiling unique alternative splicing responses to low temperature in Zoysia japonica through ZjRTD1.0, a high-quality reference transcript dataset.

Inadequate reference databases in RNA-seq analysis can hinder data utilization and interpretation. In this study, we have successfully constructed a high-quality reference transcript dataset, ZjRTD1.0, for Zoysia japonica, a widely-used turfgrass with exceptional tolerance to various abiotic stress, including low temperatures and salinity. This dataset comprises 113,089 transcripts from 57,143 genes. BUSCO analysis demonstrates exceptional completeness (92.4%) in ZjRTD1.0, with reduced proportions of fragmented (3.3%) and missing (4.3%) orthologs compared to prior datasets. ZjRTD1.0 enables more precise analyses, including transcript quantification and alternative splicing assessments using public datasets, which identified a substantial number of differentially expressed transcripts (DETs) and differential alternative splicing (DAS) events, leading to several novel findings on Z. japonica's responses to abiotic stresses. First, spliceosome gene expression influenced alternative splicing significantly under abiotic stress, with a greater impact observed during low-temperature stress. Then, a significant positive correlation was found between the number of differentially expressed genes (DEGs) encoding protein kinases and the frequency of DAS events, suggesting the role of protein phosphorylation in regulating alternative splicing. Additionally, our results suggest possible involvement of serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs) in generating inclusion/exclusion isoforms under low-temperature stress. Furthermore, our investigation revealed a significantly enhanced overlap between DEGs and differentially alternatively spliced genes (DASGs) in response to low-temperature stress, suggesting a unique co-regulatory mechanism governing transcription and splicing in the context of low-temperature response. In conclusion, we have proven that ZjRTD1.0 will serve as a reliable and useful resource for future transcriptomic analyses in Z. japonica.

KIS counteracts PTBP2 and regulates alternative exon usage in neurons.

Alternative RNA splicing is an essential and dynamic process in neuronal differentiation and synapse maturation, and dysregulation of this process has been associated with neurodegenerative diseases. Recent studies have revealed the importance of RNA-binding proteins in the regulation of neuronal splicing programs. However, the molecular mechanisms involved in the control of these splicing regulators are still unclear. Here, we show that KIS, a kinase upregulated in the developmental brain, imposes a genome-wide alteration in exon usage during neuronal differentiation in mice. KIS contains a protein-recognition domain common to spliceosomal components and phosphorylates PTBP2, counteracting the role of this splicing factor in exon exclusion. At the molecular level, phosphorylation of unstructured domains within PTBP2 causes its dissociation from two co-regulators, Matrin3 and hnRNPM, and hinders the RNA-binding capability of the complex. Furthermore, KIS and PTBP2 display strong and opposing functional interactions in synaptic spine emergence and maturation. Taken together, our data uncover a post-translational control of splicing regulators that link transcriptional and alternative exon usage programs in neuronal development.

RNA-binding proteins potentially regulate the alternative splicing of cell cycle-associated genes in proliferative diabetic retinopathy.

RNA-binding proteins (RBPs) contribute to the pathogenesis of proliferative diabetic retinopathy (PDR) by regulating gene expression through alternative splicing events (ASEs). However, the RBPs differentially expressed in PDR and the underlying mechanisms remain unclear. Thus, this study aimed to identify the differentially expressed genes in the neovascular membranes (NVM) and retinas of patients with PDR. The public transcriptome dataset GSE102485 was downloaded from the Gene Expression Omnibus database, and samples of PDR and normal retinas were analyzed. A mouse model of oxygen-induced retinopathy was used to confirm the results. The top 20 RBPs were screened for co-expression with alternative splicing genes (ASGs). A total of 403 RBPs were abnormally expressed in the NVM and retina samples. Functional analysis demonstrated that the ASGs were enriched in cell cycle pathways. Cell cycle-associated ASEs and an RBP-AS regulatory network, including 15 RBPs and their regulated ASGs, were extracted. Splicing factor proline/glutamine rich (SFPQ), microtubule-associated protein 1 B (MAP1B), heat-shock protein 90-alpha (HSP90AA1), microtubule-actin crosslinking factor 1 (MACF1), and CyclinH (CCNH) expression remarkably differed in the mouse model. This study provides novel insights into the RBP-AS interaction network in PDR and for developing screening and treatment options to prevent diabetic retinopathy-related blindness.

ASAPA: a bioinformatic pipeline based on Iso-Seq that identifies the links among alternative splicing, alternative transcription initiation and alternative polyadenylation.

BACKGROUND: Although the events associated with alternative splicing (AS), alternative polyadenylation (APA) and alternative transcription initiation (ATI) can be identified by many approaches based on isoform sequencing (Iso-Seq), these analyses are generally independent of each other and the links between these events are still rarely mentioned. However, an interdependency analysis can be achieved because the transcriptional start site, splice sites and polyA site could be simultaneously included in a long, full-length read from Iso-Seq. RESULTS: We create ASAPA pipeline that enables streamlined analysis for a robust detection of potential links among AS, ATI and APA using Iso-Seq data. We tested this pipeline using Arabidopsis data and found some interesting results: some adjacent introns tend to be simultaneously spliced or retained; coupling between AS and ATI or APA is limited to the initial or terminal intron; and ATI and APA are potentially linked in some special cases. CONCLUSION: Our pipeline enables streamlined analysis for a robust detection of potential links among AS, ATI and APA using Iso-Seq data, which is conducive to a better understanding of transcription landscape generation.

Studying the Effect of MBNL1 and MBNL2 Loss in Skeletal Muscle Regeneration.

Loss of function of members of the muscleblind-like (MBNL) family of RNA binding proteins has been shown to play a key role in the spliceopathy of RNA toxicity in myotonic dystrophy type 1 (DM1), the most common muscular dystrophy affecting adults and children. MBNL1 and MBNL2 are the most abundantly expressed members in skeletal muscle. A key aspect of DM1 is poor muscle regeneration and repair, leading to dystrophy. We used a BaCl2-induced damage model of muscle injury to study regeneration and effects on skeletal muscle satellite cells (MuSCs) in Mbnl1∆E3/∆E3 and Mbnl2∆E2/∆E2 knockout mice. Similar experiments have previously shown deleterious effects on these parameters in mouse models of RNA toxicity. Muscle regeneration in Mbnl1 and Mbnl2 knockout mice progressed normally with no obvious deleterious effects on MuSC numbers or increased expression of markers of fibrosis. Skeletal muscles in Mbnl1∆E3/∆E3/ Mbnl2∆E2/+ mice showed increased histopathology but no deleterious reductions in MuSC numbers and only a slight increase in collagen deposition. These results suggest that factors beyond the loss of MBNL1/MBNL2 and the associated spliceopathy are likely to play a key role in the defects in skeletal muscle regeneration and deleterious effects on MuSCs that are seen in mouse models of RNA toxicity due to expanded CUG repeats.

The impact of genetically controlled splicing on exon inclusion and protein structure.

Common variants affecting mRNA splicing are typically identified though splicing quantitative trait locus (sQTL) mapping and have been shown to be enriched for GWAS signals by a similar degree to eQTLs. However, the specific splicing changes induced by these variants have been difficult to characterize, making it more complicated to analyze the effect size and direction of sQTLs, and to determine downstream splicing effects on protein structure. In this study, we catalogue sQTLs using exon percent spliced in (PSI) scores as a quantitative phenotype. PSI is an interpretable metric for identifying exon skipping events and has some advantages over other methods for quantifying splicing from short read RNA sequencing. In our set of sQTL variants, we find evidence of selective effects based on splicing effect size and effect direction, as well as exon symmetry. Additionally, we utilize AlphaFold2 to predict changes in protein structure associated with sQTLs overlapping GWAS traits, highlighting a potential new use-case for this technology for interpreting genetic effects on traits and disorders.